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We study numerically the dynamics of two-qubit gates with superconducting charge qubits. The 
exact ratio of Ej to El and the corresponding operation time are calculated in order to implement 
two-qubit gates. We investigate the effect of finite rise/fall times of pulses in realization of two-qubit 
gates, ft is found that the error in implementing two-qubit gates grows quadratically in rise/fall 
times of pulses. 
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I. INTRODUCTION 

Building a practical quantum computer with a large 
number of qubits has recently attrajeted much attention 
since the useful quantwH algorithmsoH and the quantum 
error-correction codea3'cl were developed. Basic quan.- 
tum logic gates were implemented oijuixapped ionsaQ, 
QED cavitiesQ, and liquid state NMRla'M. Solid-state 
devices have advantages of large scale integration, flexi- 
bility in the design, and easy connection to conventional 
electronic devices. Some of solid-state devices-proposed 
are as follows: electron spins ixL-quantum dotaiJ, nuclear 
spins of donOT-|SutoiH|an-sU^ and ultrasmall Joseph- 
son junctionaijoHESIiZlllj. 

Two types of superconducting qubits based on ultra- 
small Josephson junctions were proposed. One is to use 
a number of excess Cooper pairs on a superconduct- 
ing ej o o pcrj pair box, called the superconducting charge 
qubitliaOllj. And the other is to utilize a single flux 
quantum of a superconducting loop, called the supercon- 
ducting flux qubitllj. Here, we focus on superconducting 
charge qubits. In theoretical aspects, the measurement 
of superconducting charge states using a single electron 
transistor and decoherence. djj e, to coupling of qubits with 
environments were studiedllZllla. Also the quantum leak- 
age of superconducting charge qubits was pointed outlia. 
In experiments, Nakamura et alS3 demonstrated the co- 
herent oscillations of Cooper pairs on a superconducting 
Cooper-pair box. This corresponds to a rotation of a sin- 
gle qubit about x axis. Thus the realization of two-qubit 
gates is an important step for making a superconducting 
quantum computer with a medium size. 

In this paper, we investigate how the two-qubit gate 
could be implemented on superconducting charge qubits. 
The pulse sequence and operation times, and the simula- 
tion of two-qubit gates are reported. We examine errors 
due to finite rise/fall times of pulses in the realization of 
two-qubit gates. 

Our paper is organized as follows. In Sec. II we intro- 
duce the Hamiltonians of superconducting charge qubits 
and of ideal qubits. A numerical method which allows 
us to simulate a quantum computer is presented. In 
Sec. Ill we explicitly show the pulse sequence necessary to 



qubits. Also numerical studies of two-qubit gates are pre- 
sented. In Sec. 1^ we analyze errors due to finite rise/fall 
times of pulses on superconducting charge qubits. Fi- 
nally, in Sec. ^ we summarize the results. 



II. HAMILTONIAN OF SUPERCONDUCTING 
CHARGE QUBITS 

In this paper we consider two qubit sjffitems: ideal 
qubits and superconducting charge qubitaljH. By com- 
paring two qubit systems, one can discover differences 
and similarities between them, which will be helpful to 
improve the design of superconducting charge qubits. 

First consider an ideal model for a quantum computer. 
The Hamiltonian of the system of N ideal qubits reads 



N 



H: 



ideal 



i=l 



implement two-qubit gates with superconducting charge 



where cr'^'-' = {ax\ <Jy\ ai^'^ ) are Pauli matrices for the 
i-th qubit. Here the Zeeman coupling terms Bi(t) and 
the inter qubit couplings Jij (t) can be turned on and off 
between zero and finite values in a controlled way. 

Let ipL.CQ]isider a system of N superconducting charge 
qubits Ejll3a. Each qubit consists of a single-Cooper- 
pair box with two ultrasmall Josephson junctions of ca- 
pacitance Cj forming a DC-SQUID ring and a gate elec- 
trode with capacitance Cg. The dynamics of a qubit is 
characterized by relevant energy scales; superconducting 
gap A, charging energy Ec = e^/2(Cg + 2Cj), Joseph- 
son coupling energy E'j, and thermal fluctuation fc^T. 
Assume that the system is in the regime of Eq <C A and 
ksT ^ Ec in order to suppress quasi-particle tunneling 
or excitation. Also suppose that the system operates un- 
der the conditions Ej ^ Ec and CgVg/{2e) ^ 1. Then 
only two charge states {|0) , |1)}, no excess Cooper pair 
and one excess Cooper pair on the box, play a role and 
represent the ottbitj—^yiew ways of coupling charge qubits 
were proposedoEjilj'ta. Here we consider the coupling 
between charge qubits via the LC resonant circuit, where 



N charge qubits are connected in parallel to a common 
inductor with inductance L. The Hamiltonian of the sys- 
tem of superconducting charge qubits is given by 

1 ^ 



i=l 

EE.TiE, 
... El 



(2) 



where Eci = 4:Ec{CgiVgi/e — 1) is turn on and off 
by applying the gate voltage Vgi on the z-th gate elec- 
trode. Here Eji{^xi) = ^E'j cos{TT^xi/^o) is the effec- 
tive Josephson energy of the z-th qubit, and controlled 
by the external flux ^xi- = h/2e is a flux quantum. 



Here is -El = 



with 



C: 



-L V Cqb ' 

The Hamiltonian of Eq. is similar to the ideal 
one of Eq. (|I]) except one fact. The system of super- 
conducting charge qubits has only two kinds of inde- 
pendently controllable parameters, {Eci^ Eji}, whereas 
the ideal model has four kinds of controllable parame- 
ters {B.j:i, Byi, Bzi, Jij}- Coupling between two super- 
conducting charge qubits could be realized by turning on 
both Eji and Ejj. This causes somewhat disadvantages 
in operating two-qubit gates, which will be discussed in 
Sec. fvl 



The time-evolution of N qubits is governed by a time- 
dependent Schrodinger equation 



ih- 1^) - H,{t) m , 



(3) 



where {ip) = X]m=o ^m{t) I™) with L = 2^ — 1 and a™ S 
C. Its computational basis is represented by a tensor 
product of individual qubits, \in) — |(7i)(8i|g2)'X'- • -(i^lqw), 
where integer m = 2^~^qi + 2^~'^q2 + • • • -l- 2*'(7jv and 
e {0, 1}. Multiplying both sides of Eq. (|) by one 
gets a set of N coupled first-order ordinary differential 
equations for the function a„(<:) 



Etnm (t) ar. 



(4) 



where -ff„m(i) = E[q{t) \m). We developed a program 
which allows us to solve the above initial value problem 
and to simulate a quantum computer. Our code is based 
on the Runge-Kutta methodn3. 

A quantum logic gate can be realized by controlling 
the time-evolution operator which acts on selected qubits 
for a fixed period of time. A way to control the time- 
evolution operator is to turn on or off each terms in 
the Hamiltonians of Eq. (jl|) or (H). For an example, a 
pulse P{t), which is applied as Bai{t) — B^^P{t) with 
a = x,y,z, makes it possible to change a term of the 
Hamiltonian from zero to finite value B^^ for a finite 
time r. In this paper, we consider a rectangular pulse 
tb — ta and a unit height, modeled 
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FIG. 1: Rectangular pulses with finite rise/fall time 2e are 
plotted as a function of time. Here time is normalized in the 
unit of nh/Bx- The role of idle time between two successive 
pulses is to prevent the tails of pulses from overlapping each 
others. 



by a kink and anti-kink pair 



tanh 



(5) 



where the rise (or fall) time is about 2e. That means 
Prccita + e) ^ tanh(2) « 0.964 for r > 2e. The smaller 
e one takes, the sharper rectangular pulse one obtains. 
Fig. 1^ illustrates a sequence of two rectangular pulses. 



III. CNOT GATE WITH SUPERCONDUCTING 
CHARGE QUBITS 

In this section, we explicitly show how the controUed- 
not (CNOT) gate with superconducting charge qubits 
could be implemented by applying a sequence of pulses. 
What makes the CNOT gate so important is that the 
CNOT gate (or any nontrivial two-qubit gate) and 
single-qubit gates form a universal set of logic opera- 
tion. Also a general two-qubit controlled-?/ gate can 
be built up of two CNOT gates and three single-qubit 
gates. The CNOT gate acting on two qubits i and 
J is represented by the unitary matrix in the basis of 

{|o.o,),|oa,),|i.o,),|ia,)} 



^CNOT 



/ 1 
10 
1 

Vo 1 



(6) 



Piecit) with width r 



where i and j are indices for a control bit and a target 
bit, respectively. 



A. CNOT gate with ideal qubits 

First consider the implementation of the CNOT gate 
for an ideal system whose Hamiltonian is given by 
Eq. (|^). The primary two-qubit gate acting on the i- 
th and j-th qubits could be implemented by turning 
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FIG. 2: Time-dependence of probabilities ja^jj^'s as a function 



of the normalized time Bxit/irh on the action of the CNOT 

of rectangular 



gate Uqnqx in an ideal model. The sequence of rectangular 



pulses is given in the below part. Here is 
with i = 1, 2. 



10 

COS 7 I sin 7 

i sin 7 cos 7 

1 



(7) 



where 7 = Jijt/h. Then, the CNOT gate for an ideal 
qubit system can be realized by combining single-qubit 
gates and the primary two-qubit gate C/j^ 



CNOT — ^ ^ 



(8) 



where is a Hadamard gate acting on the i-th qubit. 

Fig. 1^ depicts the time-evolution of qubits under the 
action of CNOT gate for the ideal model. The lower 
part of Fig. (H) illustrates the sequence of rectangular 
pulses which generates the CNOT gate of Eq. (||) . Here 
is taken j-^in) = |3) = |11), i.e., janp = 1, as an initial 
state (denoted by o). After implementing the CNOT 
gate, one obtains |aioP = 1, corresponding to state |10), 
labeled by In Fig. p|, it takes a long time to implement 
the single-qubit gate e~''~"'^ because Bai^s are taken to 
be always positive or zero. 



B. CNOT gate with superconducting charge qubits 

Let us consider the implementation of the CNOT gate 
with superconducting charge qubits. The coupling be- 
tween the i-th and j-th qubits can be switched on by 
turning on Josephson couplings, Eji and Ejj, and by 



turning off charging energy terms, E, 



a 



E, 



0. 



Then, the interaction Hamiltonian between two charge 



qubits, for an example qubits 1 and 2, can be written as 

Ej 



(2) 



SintaW^^) , (9) 



where E-mt = Ej/El with an assumption of Eji — 
Ej2 = Ej. For a moment, suppose that switching Eji'a 
on and off could be done instantaneously. As will be dis- 
cussed later, a finite time in switching on EjiS gives rise 
to an error in implementing two-qubit gates. 

The basic two-qubit gate is given by the time-evolution 
operator under the Hamiltonian of Eq. for a finite 
time r 



C^ph(r) 



(10) 



Let us transform the above time-evolution operator under 
the unitary operator Ry = exp[— z^(ct1,^'' + o'y'^'')] as 



C^ph = ^1 Uph Ry 



on the coupling Jij (t) . It can be written in the basis reads 

{|o,Oj),|oaj) |i^o,),|ia,)} as 



where the Hamiltonian transformed, H^^ 



(11) 

RtjHphRy 



^.W-^cxf -i?i„t.eV(^) (12a) 



2 2 

/a 0-1 



10 
10 
V-1 -ay 



(12b) 



Here a = Ej/Ei^t = El/Ej. Note that the Hamiltonian 
of Eq. (^2|) does not mix the subspace spanned by |00) 
and |11) with that by |01) and |10). Then the time- 
evolution operator of the transformed Hamiltonian H^^ 
becomes 



/cos 

V - 






irir sin ( 





cos 4 
i sine 






i sin ( 
cos 4 




—iTix sm t) 



cos 6 — Hz sin ( 



(13) 



where the time evolution corresponds to the rotation 
about the y-axis by angle = Ei-atT/fi in the sub- 
space spanned by |01) and |10). Also in the subspace 
spanned by |00) and |11), the time evolution gives rise 
to the rotation about the axis (na,,0,rtz) by angle 9 = 
\/l + Ei-a±T /%, where 



El 



1 Ej 



(14a) 
(14b) 



By choosing the appropriate values El, Ej, and r, one 
can control the rotation angles 4> and 6. Let us consider 
special angles 4> = f (2to — 1) with m = 1, 2, • • • , and 9 = 
niT with n = 0, ±1, ±2, • • • . Note that the case of even n 



4 



is same to that of odd n up to the global phase e*'^. These 
angles can be obtained by taking the evolution-time r = 
j{2m — 1)-^^, which is given by the relation of angle cj). 

Then one has 6* = f (2m-l)Vl + a2 = ^y/l + El/E'j = 
nTT. Thus we obtain the relation between E]^ and Ej 



(15) 



where it could be accomplished by tunning Ej or E^. 
Note that one can not take n = 0. Also the evolution- 
time is given by 




ae] 



(2m- 1) 



(16) 



We have the basic two-qubit gate Uph = Ry U!^^ 
written by 



0*0 







V 



- e 




1 



-i(f) 



1 - e**\ 









1 



, (17) 



where (j) = f (2r7i — 1) with m — 1, 2, • • • . By combin- 
ing Uph with single-qubit gates, the controlled phase flip 
gate Uqpp, operating on the i-th and j-th qubits, can be 
realized by 



CPF 



where is the value given above. By using the controlled 
phase flip gate U^p-p and the Hadamard gate Hj on qubit 
j, the CNOT gate can be implemented by 

(19) 



^CNOT ^ ^j"-^CPF^J 



Taking m 
^Ej. 



1 and 



3, one gets El 



giisi— value satisfies the physical condition 
,i8[gi| rpj^g typical time-scale for operating 
the basic two-qubit gate is given by r = -\/143, which 
is much longer than the operation time of the single-qubit 
gate, Top = fi-fEj or h/Eci- Fig. || shows time-evolution 
of qubits when \^pin) = |11) is taken as an input (denoted 
by o). After operating the CNOT gate on two qubits, 
one gets the output IV'out) = |10), which is labeled by x. 



IV. ERRORS DUE TO FINITE RISE/FALL 
TIMES OF PULSES 



Up to now, it was assumed that the effective Joseph- 
son coupling Eji{t) could be switched on (or off) in- 
stantaneously. This means a pulse applied is a per- 
fect rectangular one and has no rise/fall times. How- 
ever, in reality it takes finite times to turn on (or off) 
pulses fully. While Eji(t)^s being switched on (or off), 
the Hamiltonian of Eq. (O) Hp^ becomes time-dependent 
and the -ffph's at different times do not commute, i.e.. 




FIG. 3; Time-evolution of probabilities of qubits as a func- 
tion of time on action of the CNOT gate for superconduct- 
ing charge qubits. Time is normalized in the unit of ttTi/Ej. 
Eci/Ej = 2 and El ~ \/143>Ej are taken. The sequence of 
rectangular pulses is depicted in the below. 



[ifph(i), i?ph(i')] 7^ for t ^ t' . This causes an error in 
the realization of two-qubit gates with superconducting 
charge qubits. If Haniiltonian's at different times com- 
mute, the shape of the pulse applied is of no importance. 
The ideal qubits and single-qubit rotations of supercon- 
ducting charge qubits are_free from this problem. 

The Magnus expansioro provides the means of repre- 
senting the time-evolution operator of a time-dependent 
Hamiltonian H{t) as 

U{t) = Te-T^ fo ^^'f^(^) = e-i"^ , (20a) 

where T denotes the time-ordering operator. Here the 

average Hamiltonian is given 

where 



1 r 



dt,H{h), 



(20b) 



2Th 



Jo 



t2 



^''^ ^ TTT dh[H{t2),H{h)]. (20c) 



If the rise/fall time is short, then the approximation to 
the first term H^^^^ is good. Otherwise the higher order 
terms become significant. If the pulse for the rise time is 
modeled by P{t) = t/2e, the first term H'-^\ which arises 
when the basic two-qubit gate C/ph is implemented, can 
be written by 



5(1) « ^(1)^(2) +^(1)^(2) 



(21) 



The error due to the finite jJse/fall times of pulses is 
quantified by the gate fidelitycJ 

F = (T^inlf/VoutC/IV'in) = KV'outlV'^ut)!' : (22) 

where U is the unitary operator corresponding to the 
ideal gate when perfect rectangular pulses are applied. 



5 



On the other hand, the unitary operator generated by 
pulses with finite rise/fall times transforms the input 
state iV'in) into the density operator of the imperfect out- 
put state Pout- In our case the gate fidehty is nothing 
but the square of the overlap between the perfect out- 
put state IV'out) = U \ipin) and the imperfect output state 
Pout = IV'out) (V'outI- Assuming that Eq. ^ is small, it 
is straightforward to calculate the gate fidelity 



(23) 



where (Ajy^) = {ipin\r]'^ |V'in) - {{ipin\ V li^in))^ is the dis- 
persion of 77 = H^^'^T/he. For small rise/fall time, error 
grows quadratically in rise/fall time 2e. 

(a) Ideal qubits 
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Fig. ^ shows the numerical study of the error due to fi- 
nite rise/fall times in implementing the CNOT gate with 
superconducting charge qubits. For ideal qubits, if fi- 
nite idle times between successive pulses are applied, the 
correct CNOT gate is realized. However, for supercon- 
ducting charge qubits, although a finite idle time reduces 
the error, there still exits the error due to finite rise/fall 
times. As shown in Fig. ^ (b), the probability of getting 
the correct state reduces quadratically in rise/fall time 2e 
in agreement with the theoretical prediction of Eq. (122 



In Nakamura et a/.'s experimentnJ, the rise/fall times 
of the pulse was about 30 — 40 ps at the top of the cryo- 
stat. If the effective Josephson coupling is Ej w 50 fJ.eV, 
then the timescale of single-qubit x rotation is about 
h/Ej w 1 ps. This means that the rise/fall times of 
the pulses should be less than 1 ps in order for two-qubit 
gates to be implemented correctly. It should be noted 
that this type of errors is caused by the coupling scheme 
between two qubits. Thus this problems can be solved 
by improving the design of devices, i.e., by introducing 
new ways of two-qubit couplings. 



SUMMARY 
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We have studied the dynamics of qubits on the ac- 
tion of CNOT gates for ideal and superconducting charge 
qubits. We have explicitly shown how the CNOT gate 
could be implemented for superconducting charge qubits. 
It is found that the error in implementing two-qubit gates 
with superconducting charge qubits grows quadratically 
in finite rise/fall times. Thus it is necessary to keep the 
rise/fall times small or to find new ways of coupling two 
qubits other than the coupling scheme via the common 
inductor. 



FIG. 4: Probability of success on the action of the CNOT 
gate as a function of the ratio of the rise/fall time to the 
operation time, e/rop: (a) ideal qubits (b) superconducting 
charge qubits. Here the time scale of operation is Top — h/B^ 
for the ideal case and Top — h/Eci for the superconducting 
case. Idle times prevent the tails of pulses from overlapping. 
The solid line is the fitting function 1 — 0.027(e/rop)^ which 
shows the quadratic growth of error in rise/fall times. 
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